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We investigate the effects of several Abelian Aharonov- 
Bohm fluxes cj> on the Euclidean Dirac spectrum of light 
quarks in QCD with two colors. A quantitative change in 
the quark return probability is caused by the fluxes, resulting 
into a change of the spectral correlations. These changes are 
controlled by a universal function of <tl0 2 where <tl is the per- 
tinent Ohmic conductance. The quark return probability is 
sensitive to Abelian flux-disorder but not to Z% flux-disorder 
in the ergodic and diffusive regime, and may be used as a 
probe for the nature of the confining fields in the QCD vac- 
uum. 

PACS numbers: 11.30. Rd, 12.38.Aw, 64.60.Cn 



1. In a finite Euclidean volume V light quarks exhibit 
properties analogous to that of electrons in small metallic 
grains Jj|]. The finite volume allows for a discrimination 
of length scales with the emergence of a diffusive regime 
characterized by the diffusion constant D w 0.22 fm [Q. 
Although QCD confines, the light quarks diffuse in four 
dimensions as a result of the spontaneous breaking of 
chiral symmetry |Q,0|. 

In two-color and two-flavor QCD the quarks are in 
the pseudoreal representation of the flavor group 
and the spontaneous breaking of chiral symmetry is ac- 
companied by the occurrence of five Goldstone modes: 3 
pions and a baryon and an antibaryon. The five Gold- 
stone modes decay weakly, with a mass that satisfies the 
Gell-Mann-Oakes-Renner (GOR) relation ||. The latter 
mass vanishes in the chiral limit. 

The nature of the Goldstone modes and the character- 
istics of the light quark spectrum are intimately related 
in the general context of chiral disorder 0. In particular, 
two-color QCD is characterized by diffuson (diagonal, pi- 
onic) and cooperon (interference, baryonic) contributions 
in the semi-classical analysis. The cooperons correspond 
to weak-localization (coherent back-scattering) contribu- 
tions to the quark return probability They are 
soft because the baryons in two-color QCD are diquarks 
and soft. QCD with three colors admit only diffusons as 
the diquarks are expected to be heavy. As well known 
in disordered metals |(J, the weak-localization contribu- 
tion is altered by the presence of an external parameter 
that breaks time-reversal symmetry. An example is an 
Aharonov-Bohm flux. 

In this letter we consider the effects of several Abelian 



Aharonov-Bohm fluxes on the light quarks in a finite Eu- 
clidean volume with two-color QCD. In many ways, our 
analysis in two-color QCD will parallel recent analyses 
of electrons in disordered metals ||. We discuss these 
effects on the light quark return probability, and analyze 
their importance on the spectral correlations. The fluxes 
cause a periodic interpolation in the spectral statistics be- 
tween orthogonal and unitary ensembles. We show that 
the quark return probability is sensitive to Abelian flux- 
disorder but not to Z2 disorder, and discuss the relevance 
of this result for the nature of the QCD vacuum. 

2. In a finite Euclidean volume V — L A pierced 
by Abelian (essentially electromagnetic) fluxes M = 
(01, 02, 03, 04) , the light quarks in the background gluon 
field A satisfy the following Dirac equation 



if[A}q k = \ k [A,(f>]q k . 
subject to the boundary condition 

q k (x + L„) = -e i2 ^ q k (x) 



(1) 



(2) 



The 0's are given in units of an electromagnetic flux 
quantum 2h/e |?J set to 1 for convenience. Through (|^) 
the quark spectrum (|l|) depends explicitly on 0. We note 
that the Abelian fluxes through (|2|) are simple holonomies 
along each of the four directions of the four- volume: 
e i2Trx-4>/L^ rp ne e ig enva i ue problem (|l[|3) is equivalent to 
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if[A] + — 7 • 



I qk = A fe [A,0] q k 



(3) 



with anti-periodic boundary conditions. 

The probability p(t, 0) for a light quark to start at a;(0) 
in V and return back to the same position x{t) after a 
proper time duration t, is 

P(*,0) = ^(l(^(0)|e 4 ™ +ml)|t| |x(0))| 2 ) A . (4) 

The averaging in (^) is over all gluon configurations using 
the unquenched two-color QCD measure. The normaliza- 
tion in (|j) is per state, where N is the total number of 
quark states in the four-volume V. Equation (|4|) may be 
written in terms of the standard Euclidean propagators 
for the quark field, 
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x^Tr (S(x,y; zi, 4>)S^ (x,y; z 2 , <f>)) ^ (5) 

generalizing to 4> t nc results in jlj. Here Z\,i — 
m — i\\ 2, and 



S{x,y;z,<f>) = (x\- 



1 



■ iy[A] + ■ 
subject to the boundary condition 

S(x + L ll ,y;z,<f>) = -e*"*'> S(x,y;z, 



(0) 



(7) 



Setting Ai,2 = A ± A/2 and neglecting the effects of A 
in the averaging in (|5j) (zero virtuality), the correlation 
function in (|^) relates (in general) to an analytically con- 
tinued pseudoscalar correlation function, as the eigenval- 
ues qk and 75 qk are pair-degenerate (chiral) jjj]. 

For two-color QCD there is an extra symmetry 
that makes the pseudoscalar correlation function degen- 
erate with certain diquark correlation functions in the 
flux-free case. Indeed, for <fi = and K = —T 2 CK the 
eigenvalues qk and K qk are also pair-degenerate. Here C 
is the charge-conjugation matrix, T 2 the color matrix and 
K the (right-left) complex-conjugation. For instance, for 
two flavors the three pions are degenerate with the spin- 
zero baryon and antibaryon. 

In general, we may rewrite (0) in the form 



p(t, 



EV 



2ttN y- 



lim 



dX 
2^ 



e~ lX ^C G {x,y-z^) (8) 



with z = m — i\/2 and E = J dA. For <j> = 0, the above 
symmetries allow us to write (two flavors) 



C G (x,y;z,0) = 

+ {S{x,y;z,0)ij 5 T 2 S(y,x;z,0)ij 5 T 2 / , 



(9) 



The flavor matrix t 2 was introduced to retain the con- 
nected parts in the correlator version || . The first contri- 
bution is pionic (diffuson), while the second one is bary- 
onic (cooperon). Both are bosonic and identical. The de- 
composition (M is commensurate with the semi-classical 
description (long paths) where the time-reversed trajec- 
tories are retained Jl]||. For 0^0, the first symmetry 
(chiral) is retained while the second one is upset. Indeed, 
now the contributions in (||) arc no longer the same as 
the fluxes add in the cooperon, but cancel in the diffuson. 
For long paths, we have 



C G (x,y;m,< 



,iQ-{x-y)_ 



1 



,iQ-(x-y) 



F 2 Q 2 + Til 2 , 

s 2 1 

F 2 Q2 + m 2 



(10) 



with = n^2n/L and = (n^ + 2(j) ^2-k / L . Here, 
£ = \{qq)\ and F is the weak decay constant for the 
Goldstone modes of mass mg, each of which are flux 
independent to leading order. 

Using the GOR relation F 2 m 2 G — mS and the analyt- 
ical continuation m — * m — iX/2, we find 



C G {x,y]z,4>) fa + 7^7 E 



2E 



-i\+2m+DQ 2 



1 V p^-O-jO — 

2V^ 



Q 



-i\+2m+DQ 2 



(11) 



with the diffusion constant D = 2F 2 /£. Inserting jll| ) 
into (|), and noting that E/A = N and p = 1/AV, with 
£ = 7T/9, we conclude after a contour integration that 



Q 



,-DQ 2 \t 



(12) 



The cooperon contribution is periodic in the flux <f> with 
periodicity 4> = 0, ±1/2, ±1, and reflects on the fact 
that the soft part of the spectrum does not discriminate 
between bosonic or fermionic boundary conditions in the 
flux-free case. The cooperon contribution may be rewrit- 
ten using Poisson's rcsummation formula as 



pc(t, 4>) 



V 



2(4nDt) 2 ^ 



Y^ e -2m|t|- 



cos (4^^) (13) 



with integer Vs. This result is in agreement with the one 
derived by Montambaux in disordered metals in lower 
dimensions. The flux-accumulation in the cooperon part 
implies changes in the spectral correlations of the light 
quarks as we now show. 

3. To describe the spectral correlations associated with 
(|l|) we will use semi-classical arguments for the two-point 
correlation function R(s, cj>) of the density of eigenval- 
ues |l]||[§. Its spectral form factor K(t, <p) is defined as 



R(s,< 



dti 



sAt 



K(t,< 



(14) 



where A = 1/pV is a typical quantum spacing at zero vir- 
tuality and in the absence of a flux. For diffusive quarks 
in two-color QCD, K (t, (f>) relates to the return probabil- 
ity through 



K(t,<f>) 



A 2 jt] 

(27T) 2 



P (t, </>) 



(15) 



which is periodic with periodicity </> = 0, ±1/2, ±1, .... 

For cj) <C 1, the flux periodicity may be neglected, and 
for times larger than the ergodic time r crg = L 2 / D but 
smaller than the Heisenberg time tjj = 1/A, the dom- 
inant contribution to the return probability stems from 
the zero modes in (0). Hence 
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with 



16tt 2 
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16"7T <7£ 



(16) 



(17) 



Here ol, = D/(AL 2 ) = 2(FL) 2 /tt is the dimensionless 
Ohmic conductance In this limit, both the return 
probability and the spectral form factor are only a func- 
tion of the combination ol</> 2 . The result (0) could be 
qualitatively derived by noting that if 2 x 2n(f> is the typ- 
ical accumulated flux per winding path |]l0f in the quark 
return probability, then I windings generate an accumu- 
lation $ = Since the paths are diffusive with an 
arc length s = LI ~ VTtt, then on the average 



Hence 



<$ 2 ) = <( 47 r/0) 2 )^(4^) 2 ^ 



( e »\ „ e -<**>/2 = e -^l*/tel 



(18) 



(19) 



which is the dephasing appearing in ( [Tq ) between the 
diagonal and interference terms. An average 2tt accumu- 
lation in $ results from quark trajectories with a long 
proper time t ~ 1/E C where E c = D/L 2 is the Thouless 
energy [|lj. 

The spectral rigidity around zero virtuality (and also 
in bulk) follows from j@|J 



E 2 (7V,0) 



-N 



ds(N- \s\)R(s, 



— N 



with N = E/A > 1. In particular 
1 



E 2 (AT,<£) 



2tt 2 



rln 



1 



TV 2 

^2" 



iV 2 



(20) 



(21) 



with a = 2m/ A and a = a + F 2 . For a ^> 1 the spectral 
rigidity is reduced by about a factor of 2. It is brack- 
eted by the orthogonal (no flux) and unitary (finite flux) 
results. 

For large fluxes, the non-zero momentum contribu- 
tions have to be retained to enforce the proper = 
0, ±1/2, ±1, ... flux-periodicity. For simplicity, consider 



the case 



(0, 0, 0, 4>) with only one-flux retained. For 



long proper times, the zero modes along the 1, 2, 3 direc- 
tions contribute only, giving 



p(t,<t>) = \e-* nM x 



-4tt 2 (n+2<p) 2 a L \t/t 



") 



(22) 



which is the analogue of a diffusion in d=l. For t < r< 



org 



L /D the diffusive paths are short and do not accumulate 
enough flux. The spectral rigidity in this case is just 



E 2 (JVT,0) 



2tt 2 



TV 2 



N 2 



(23) 



with a n — a + 4ir 2 aLn- 2 and a n = a 
For N, er £ > a, (53) simplifies to 



- An 2 a L (n + 24>) 2 . 



S 2 (AT,0) = £ 2 (JV,0) 



1 



rln 1 



4— sin 2 (2^)1 (24) 
a / 



in agreement with a result derived by Montambaux || 
in the context of disordered metals. We note that for 
a small flux (|24|) is in agreement with (pl|). Interest- 
ingly enough, the conductance cr^ of the chiral vacuum 
is directly accessible from the spectral rigidity through 
( p4| ) providing for a direct measurement of this impor- 
tant quantity in disordered QCD. (53) may be assessed 
using current lattice QCD simulations. 

4. Since the quark return probability and the spec- 
tral rigidity are sensitive to flux- variations in a finite Eu- 
clidean volume, they could be used to probe the flux- 
content of the two-color QCD vacuum - in particular 
monopole-antimonopole rich vacua which are likely to 
generate flux-disordered vacua. Although our analysis 
so far has focused on Abelian electromagnetic Aharonov- 
Bohm fluxes, we speculate that in the maximally pro- 
jected gauge, colored magnetic monopoles and anti- 
monopoles generate colored and Abelian fluxes that act 
randomly on the light quarks in the vacuum. 

If we consider an Abelian flux-disordered vacuum char- 
acterized by a Gaussian distributed flux with a mean 



(25) 



then the quark return probability can be easily estimated 
from ( ]l2| ) using the Poisson form (13) for the cooperon 
part. If we split the quark return probability p = pd+PCi 
then the diffuson part po is flux- insensitive 



«p D (t,0» = 



V 



,-2m\t\ 



2{A-KDt) 2 

while the cooperon part pc is flux-sensitive 
V 



«Pc(t, 0))) 



2(4vrL>t) ; 



E 

[2] 



-2ro|t| 



(26) 



• (27) 



We note that the periodicity in <p — 0, ±1/2, ±1, ... of 
the quark return probability implies that the latter is 
likely insensitive to a Z 2 flux-disordered vacuum. If these 
effects extend to an Abelian flux-rich vacuum, a simple 
way to detect them is to measure the relative ratio of the 
quark return probabilities for N c — 2 and N c = 3. A 
flux sensitivity implies a t-dependent ratio close to 1/2 
(as opposed to 1), assuming that the vacuum structure 
does not change appreciably from two to three colors. 

5. Finally, since the transition from an orthogonal to 
unitary ensemble sets in the ergodic regime with t > r org , 
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the ensuing spectral statistics are amenable to an analysis 
in O-dimension (matrix model). The simplest realization 
that is chiral and embodies the essentials of a single flux 

is 



im A + itpH/y/N 

A + icpB/y/N im 



(28) 



where A and B are real N x N symmetric and antisym- 
metric random matrices respectively, with a fixed vari- 
ance S 2 /N . For ip = and m = the matrices are chiral 
orthogonal, while for <p ^ they are in general chiral 
unitary, hence the interpolation. The 1/VN in the off- 
diagonal matrix elements follows from second order per- 
turbation theory since the typical level shift induced by 
such terms is 



AE 



E 



E n — Eq 



A,B 



t 

N 



(29) 



where the \I/ n 's are generic eigenfunctions of (28) with 
eigenvalues E n for B — 0. In the mesoscopic (ergodic) 
limit, the sum is dominated by the modes where the en- 
ergy denominator is of order A ~ 1/N . If the states are 
delocalized then the eigenfunctions are of order l/y/N. 
Since the sum over I, J is of order iV 2 , and the aver- 
aging over B is of order E 2 /7V, the second order shift 
is AE ~ 1/7V and comparable to the mesoscopic level 
spacing A ~ 1/N. This result is consistent with the ob- 
servation made using non-chiral matrices |L4| . We recall 
that in d ^ dimension, the corresponding terms are 
down by 1/ y/V as sugg ested by (§). 

6. We have shown that in two-color QCD, Abelian 
Aharonov-Bohm fluxes cause changes in the quark return 
probability with important consequences on the spectral 
form factor and the spectral statistics. The change is 
controlled by a universal function of n 2 <jL<t> 2 , where n is 
the number of applied fluxes. This analysis offers a sim- 
ple measurement of the Ohmic conductance <jl- When 
combined with the value of the chiral condensate E, this 
leads to an independent assessment of the pion weak- 
decay constant F and a simple vindication of the GOR 
relation in two-color QCD. 

The sensitivity of the quark return probability to 
Abelian fluxes for two-color QCD raises the interesting 
possibility of addressing the nature of the confining fields 
in the QCD vacuum. Indeed, the quark return probabil- 
ity for two colors can discriminate between a Zi flux- 
disordered phase and an Abelian flux-disordered one, as 
it is blind to the former. Also, if the flux-disordering 
is generic in two and three colors QCD, quantitative 
changes in the proper time behavior of the pertinent 
quark return probabilities and spectral rigidities are ex- 
pected. The present results can be tested by using lattice 
QCD simulations, or continuum models of the QCD vac- 
uum lEHI. 
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